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Recent results on second order quantum corrections to scmileptonic decays of b 
quarks and positronium spectroscopy are reviewed. Similarities between calcula- 
tions in these seemingly different problems are demonstrated. Technical aspects of 
solving systems of recurrence relations, and application of dimensional regulariza- 
tion to bound states are explained. 



1 Introduction 

Precision measurements play a significant role in contemporary particle physics, 
spanning the energy scales starting with the Z-boson pole at the high end, 
down to scales of the order of the b and c quark masses, and further down 
to hydrogen, positronium, and muonium spectroscopy, and anomalous mag- 
netic moments of leptons. The role of the precision measurements reflects the 
successes achieved by the Standard Model (SM) of electroweak and strong in- 
teractions. It has become clear in recent years that finding any signal of "New 
Physics," beyond the SM expectations, is a very difficult task. 

Searching for hints about a more fundamental theory, which would replace 
the SM, but lacking devices necessary to make the next step up in energy, one 
resorts to a more precise analysis of the already available data and tries to 
increase the accuracy of the experiments at achievable energies. As a conse- 
quence, more and more accurate theoretical predictions are required practically 
at every point of the energy scale up to ~ 100 GeV. 

Theoretical studies are often impeded by complications inherent to bound 
states or QCD uncertainties. Consider for example the inclusive decay width 
of the i3-meson, T{B Xccv), used for the \Vcb\ determination. To calculate 
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the decay rate in a model-independent way, one first has to clarify how to 
treat the non-perturbative bound state effects. This non-trivial issue is solved 
by invoking the Heavy Quark Expansion which permits a calculation of the 
inclusive semileptonic decay width as a systematic expansion Id in powers of 
Aqcd/w. The lowest order of this expansion corresponds to the free quark 
decay. There remains a technical difficulty of performing the standard pertur- 
bative calculations. Until recently, only one- loop accuracy was achievable for 
charged particle decays. Our approach to performing two-loop calculations is 
one of the subjects of this talk. 

A complicated situation occurs also at the very low energies, where Couloni- 
bic bound states like hydrogen, muonium, or positronium are studied. The 
measurements in such systems are very precise. The theoretical predictions, 
however, are complicated by the essentially non-perturbative nature of bound 
state calculations. Fortunately, it turns out that the properties of the non- 
relativistic bound states, such as energy intervals and lifetimes, can be de- 
scribed in terms of series in powers and logarithms of the fine structure con- 
stant. Here a proper organization of the calculations is a challenging problem. 

It is common in theoretical physics that problems can be solved if a small 
parameter is available. In case of the "standard" perturbative calculations the 
small parameter is usually the coupling constant. For the bound states, the 
smallness of a coupling constant is insufhcient, because the very existence of 
the system can only be described in terms of an infinite number of Feynman 
diagrams. The situation is simplified in case of non-relativistic bound states, 
which encompass the light atoms and "leptonia" for which precision measure- 
ments are being performed. There, the relative velocity of the bound state 
components is small and can serve as an expansion parameter. 

For the bound states, expanding around a non-relativistic limit is the only 
way known at present to arrive at precise results. But even in case of the more 
conventional perturbative calculations, Feynman integrals often are intractable 
and the only option to evaluate them is to invent a small parameter. It is 
often convenient to use some kinematic parameter, so that gauge invariance 
and other symmetries of the problem are preserved. 

Construction of non-conventional expansions of Feynman graphs has re- 
cently helped to solve many outstanding problems in high precision theoretical 
studies. In what follows we are going to describe some examples of how such 
expansions can be constructed, and review our recent calculations devoted to 
the semileptonic decay width of the b quark and the positronium energy levels. 
We deliberately focus on the technical aspects of those calculations, to provide 
additional information and supplement our journal publications. 
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Maximal Recoil Line mass 



Figure 1: Status of the two-loop QCD corrections to the decay b — > c + leptons. The dashed 
line denotes the physical region for the actual c quark. Points where the full correcti«HS are 
known are circled. An analytical formula is known along the whole zero recoil line. Q'Q The 
other two pointe are found from expansions around the bpKC points of the two arrows: at 
maximal recoil Q and at the intersection with the diagonal. Q 



2 Inclusive semileptonic decay width of the b quark 

Let us start by discussing a calculation of the second order QCD corrections, 
0{a^), to the semileptonic decay width of the b quark, r{b ^ cev^. At present 
a complete calculation is not feasible. However, one can calculate corrections 
at this order to the differential decay width dF/dq^^ for special values of the 
leptonic invariant mass, q^^^ = q^. Such calculations were performed at three 
points: 

• = (jTLh — TOc)^ (zero recoil);Bi 

• = ml (intermediate recoil) ;i 

• — (maximal recoil) . @ 

These kinematical configurations are depicted in Fig. |^. Based on these 3 
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results, we performed a fit to interpolate the corrections for arbitrary . From 
that fit, a correction to the inclusive decay rate, r(6 — > cev^), was estimated. 

For kinematical configurations away from the zero recoil limit one should 
consider three different types of corrections: pure virtual corrections, single 
gluon radiation, and radiation of two gluons (at zero recoil only the first type is 
present). Our approach is to calculate these contributions separately and then 
sum them up in the final result.^ Below we describe the essential ingredients 
of those calculations. 



2.1 Virtual corrections 

As an example of the virtual correction calculation we consider the intermedi- 
ate recoil case, equating the invariant mass of the lepton system, to the 
c-quark mass. If the c quark had exactly half the h quark mass, this configu- 
ration would coincide with the zero-repj^l limit (see Fig. [^) , in which case the 



virtual corrections are known exactly. U'Q For the actual c and b masses these 
corrections can be calculated by expanding the amplitude in power series in 
{mb — 2mc)/mb, that is around the known zero recoil case. An important aspect 
of this expansion is that it commutes with the virtual momenta integrations; 
it simply amounts to Taylor expanding the integrand. 

After the expansion, all integrals we need are of the zero-recoil type. They 
can be divided up into the following categories (they are written here in a 
general form, with lu denoting the ratio of the final and the initial quark masses. 
For our intermediate recoil calculation we need only uj — 1/2): 

1. Planar with five propagators 

A(ai, a2, as, 04, as) = t:^ 
^ f d^fcid^fca 



2. Non-planar, first type, five propagators 

d^fcid^fca 



iVp(ai,a2,a3,a4,a5) = tt 



D 



l2ai 72a2 
'^l "'2 

1 



(fc2 + 2kipr (M + 2k2pujr [(fci + + 2(fci + k2)pT 



(2) 



''An alternative approach to such calculations is based on the optical theorem; to get an 
n-loop result, one calculates the imaginary part of a certain (n -|- l)-loop Feynman iirtegral. 
Such method was applied e.g. to calculate QED corrections to the muon decay rate. Q 
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3. Non-planar, second type, five propagators 



iVp(ai, 02, as, 04, as) = tt 



2 



y^2ai ^2a2 



2 

(3) 



+ 2fcip)"-' (fc2 + 2A:2pcc;)''" + fca)^ + 2cj(fci + k2)p\ 
4. Non-planar with six propagators 

d^fcid^fca 



^6(01,02,03,04,05,06) = TT 



1 



[(fci + k2)^ + 2{ki + k2)pr [(fci + + 2t<j(fci + k2)pr m + 2cuk2PT' 

(4) 

Taylor expansion of the amplitude may result in very many terms, differing 
by powers of the propagators, o^. It is crucial to find a general algorithm of 
calculating them, so that an algebraic manipulation program can be employe^ 
It is convenient to use recurrence relations derived using integration-by-parts 
to reduce an arbitrary integral to a small basis set of the so-called master 
integrals. 

Let us note at this point that solving recurrence relations is still more of an 
"art" than "science." It would be very desirable to have a better understanding 
of the mathematical structure of the systems of recurrence relational a very 
promising approach to this problem is being developed by P. Baikov.u 

Here, we would like to demonstrate our approach to solving such recurrence 
relations with the example of the planar integral Psdo^}). First, a system of 
relations is derived in the usual wayQ by equating integrals of total derivatives 
to zero. One finds 

Ml = D - 2ai - 02 - 04 + 022+(3" - 1") - 044+I" 
M2 = L» - 01 - 2o2 - 04 + oil+(3" - 2") + 044+ ('^^—^3- -2- + -5- 

M3 = D-ai-a2- 2o4 - oil+4" + 022+ f ^3" - 4" + — ) + 2a44+ 

\ UJ LO J 

M4 = -D - 202 - 03 - 05 -f 033+(l" - 2") + 055+ ((1 - w)l" - 2" + LoA-) 
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2D — 2ai — 2a2 — 2a3 — 04 — 05 — 044 1 



a; 




Me 



I? - 02 - as - 205 + 022+ ((1 - Lj)l +uj4 -5 ) 
-038+5" + 2w^a55+ 



M7 



( 



3D 
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Gi (1- - 4-) + 2a44+l- + a55+cj(l- - 2" + 3") 




3i:> 
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(3" - 5") + a44+(l" - 2" + 3") + 2cja55+3 



(5) 



The recursive algorithm proceeds in three steps: 

1. the exponents of the massive propagators (04^5) are reduced to 1; 

2. the exponent of the massless propagator containing both integration im- 
pulses (02) is reduced to 1; 

3. the exponents ai^a are reduced. 

The first step is simple. We assume that 04 and 05 are positive; otherwise 
the integral can be done by a sequence of one-loop integrations. We use the 
relation Afg. The term 21^^055+ involves an integral with 05 higher than any 
other term in this relation. Therefore we can use Mq to express P5(. . . , 05 + 1) 
by P5(. . . ,05). An analogous manipulation is done using to reduce 04. 
Eventually we obtain 04, 05 < 1. 

The next step is accomplished as follows. We first would like to eliminate 
all operators that contain 4+ and 5+ from the system (||) . There are eight such 
operators 4+, 4+1", 4+2", 4+3", and5+, 5+1", 5+2", 5+3". An important 
point is that this number is equal to the number of the recurrence relations 
and we can therefore solve the recurrence relations for these operators. 

We denote the resulting 8 equations by T4^4i_42,43,5,5i,52,53. The notations 
are such that for example T41 means an equation which expresses the operator 
4+l~ through other operators. These equations allow us to obtain new recur- 
rence relations using a4(T4i — 1"T4) = and similar equations. In particular, 
we use a relation a4(T43 — 3~T4) = to lower the value of 02 to 1. 

After this has been achieved, we eliminate all operators which contain 2+. 
At this stage, the problem is essentially two parametric (i.e. we have to move in 
the oi — 03 plane). Finally, we are able to reduce any ^5({ai}) integral to four 



integrals ^5(1,1,1,1,1), P5(0, 1, 1, 1, 1), P5(-l, 1, 0, 1, 1), and P5(0, 1, 0, 1, 1) 
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and a number of trivial integrals. The four master integrals are 




In^ 2 157r2 



2 32 
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The above approach to solving recurrence relations is quite general. One sees 
that such solution cannot be constructed in a completely universal way. How- 
ever, our approach allows to go through a number of steps, constructing a 
solution without much effort and it has proved useful in various cases. In 
particular, it works for the "hard" threshold integrals, necessary for the calcu- 
lations described in Section 

2.2 Emission of two gluons in semileptonic b decays 

For simplicity, we will discuss here the case of the maximal recoil, — 0. 
The intermediate recoil case is treated in an analogous way. A more detailed 



summarized. 

The basic observation in this part of the calculation is that the energy losses 
caused by radiation off a non-relativistic particle are small. This is evident in 
the Coulomb gauge, since the interaction of magnetic photons with charged 
particles is proportional to the velocity of the latter. Therefore, the physical 
properties of the real radiation provide a hint how an expansion should be 
organized. Considering a decay b ^ c + e + + gi + 92 in the case = {e + 
i^e)^ = 0, we intend to calculate it expanding in (5 = {mb — rnc)/mi,. For brevity 
we will refer to the e, i/ system as to the W boson. Taking into account phase 
space constraints, it is easy to establish that the spatial momentum of the c 
quark, as well as the four-momenta of the W, gi and g2 are of the order S. Then, 
a natural expansion of the propagators emerges: all propagators of the virtual 
particles can be expanded around the static limit, which essentially simplifies 
the calculation. Clearly, such an expansion will generate a large number of 
terms, but this can be handled using symbolic manipulation programs. 

The last non-trivial question is the integration of the obtained expression 
over the four particle phase space. Given that the remaining propagators are 



discussion of this issue is given 




Here the most essential features are 
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static and the W boson does not participate in the interaction, it is clear 
that the integration over two-particle phase (sub)space {W, c) can be easily 
performed. The remaining part of the calculation requires an integration over 
the three particle phase space - two gluons and the "particle" composed of W 
and c. A systematic expansion of this object is accomplished by an appropriate 
choice of variables, described in detail inE3. 



2.3 One gluon radiation and one virtual correction 

Single gluon emission in the presence of a virtual loop is, in principle, possible 
to compute exactly. However, it is very inconvenient, e.g. because of necessity 
of integrating four-point functions over a three particle phase space. 

We have adopted a different approach, based again on an expansion of 
such amplitudes around the static limit {mf, — rric). In contrast to the virtual 
corrections and double gluon emission, in the present case the expansion is non- 
trivial. Let us consider the following diagram: first a virtual gluon is exchanged 
between c and h quarks, and then the c quark emits a real gluon. An essential 
point is that before the real gluon emission the c quark is off mass shell, which 
provides an infra-red regularization for the virtual loop. Expanding the one- 
loop amplitude naively, i.e. putting the "intermediate" c-quark line loop on 
shell, we change the infrared properties of the diagram and divergences will 
result in singularities in 1/(D — 4) (we use dimensional regularization). To 
avoid these spurious singularities tlHLexpansion should be "corrected" ; in this 
case the so-called eikonal expansionEJO is used. In the context of the two 
corrections to fermion decays this expansion is described in detail in Ref. 



2.4 General remarks on semileptonic heavy fermion decays 

Let us summarize the technical experience we have gained by studying semilep- 
tonic b and t decays to second order in the strong coupling constant. 

First, we note that the problem wasjuade tractable by inventing a "small" 
parameter. For the top quark decays Ej the result is obtained if we take the 
"small" parameter close to unity. This clearly pushes the applicability of the 
method to its limits. However, even in this extreme case reasonable numerical 
estimates can be obtained. 

Second, because we consider processes with the on-shell external particles, 
one has to be careful in analyzing contributing regions of integration momenta. 
In particular, one can encounter non-analytical dependence on the expansion 
parameter. In many cases the relevant momentum regions can be found by 
considering the corresponding effective theories. 
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Third, it is possible to construct an expansion of the ampHtudes with real 
radiation of massless particles starting from the non-relativistic limit. It is 
likely that this construction can be applied beyond the semileptonic decay 
width, for example in a perturbative calculation of the cross section of such 
process as e^e~ — > fj,'^ fi~^j close to the two muon threshold, or even for higher 
energies if one takes into account sufhciently many terms of such an expansion. 



3 Bound States in QED 

Consistent description of bound states is one of the most complicated prob- 
lems in Quantum Field Theory. The problem is simplified if the bound state 
components are non-relativistic. This is the case, for example, in QCD in T or 
ip meson families, or in QED in light atoms and in leptonic bound states like 
positronium or muonium. 

The difficulty in the study of bound states is their essentially non-per- 
turbative nature, as we discussed in Section |l|. In the non-relativistic case 
the situation is simplified: to first approximation one can describe the bound 
state using a Schrodinger wave function. This description can subsequently be 
refined by means of perturbation theory, using the relative velocity v of the 
components as an expansion parameter. 

This kind of expansions can be constructed both for QED and QCD. In 
what follows, we limit ourselves to the QED where a low-energy effective the- 
ory approach is well establiahed, and is known as Non-Relativistic Quantum 
Electrodynamics (NRQED).0 



3.1 NRQED 

Very roughly speaking, the idea underlying NRQED consists in separating 
the soft scale ~ ma, which characterizes the bound state, from the hard one 
'-^ TO, and "integrating out" the latter. The interactions at the soft scale can be 
described by Quantum Mechanics, i.e. employing a non-relativistic approxima- 
tion. The physics at hard scales is governed by the relativistic QFT. However, 
when we are interested in processes which occur at small momenta and energies, 
the contribution of hard scales can be accounted for by introducing effective 
local operators and computing their coefficients in the effective Hamiltonian. 

In the lowest order of this construction we describe the system by a bound 
state wave function. For positronium it is found from the Schrodinger equation 
with the "leading order" Hamiltonian (to is the electron mass, r is the distance 
between the electron and the positron): 
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In a 3-dimensional space the ground state wave function (in momentum rep- 
resentation) is 



, , , I nam 2rr?o? mo? , , 

The difference between the full Hamiltonian and Eq. (||) can be treated as 
a perturbation. The difficulties of calculating in NRQED are caused by di- 
vergences in-the matrix elements of the subleading operators. In the original 
formulationta the ultraviolet divergences were dealt with by subtracting coun- 
terterms defined with help of diagrams evaluated on-shell and at threshold. 
Those diagrams, in turn, contain infrared divergences, which are regularized 
by introducing a small mass A for the photon. 

Since that pioneering work it has become clear that working with a massive 
photon leads to technical difficulties, particularly in evaluating the so-called 
hard corrections, arising at scales ~ m. Indeed, the hard corrections remain 
unknown (at two-loops) for such processes as para or orthopositronium decays, 
while the soft scale effects have already been evaluated. The reason for this 
is that the two-loop diagrams are difficult to evaluate when more than one 
mass scale is present. Although the only relevant scale in hard corrections for 
positronium physics is the electron mass, regularization with a photon mass 
introduces a second scale A. The hard corrections can, however, be computed 
in many configurations, if instead of introducing A one uses dimensional regu- 
larization. 

It has been shown Ej that dimensional regularization permits an exact 
separation of effects arising at various characteristic energy scales. pLIsing that 
method, the complete energy spectrum of Ps has been reproducedllj to order 
m(X^ . Further, we have computed analylically mc^ corrections to the hyperfine 
splitting (HFS) of the Ps ground state. EJ More recently,E3 we generalized that 
result to all S states and computed also their spin independent shift at 0(moP\ 
Here we would like to discuss some features of the NRQED in dimensional 
regularization, which we found useful in those calculations. 

3.2 Quantum Mechanics in d dimensions 

The first point we would like to discuss is a construction of the lowest or- 
der approximation. We consider the Schrodinger equation in a d-dimensional 
momentum space, 

Anam f d^k 0(fc) 
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The solution of this equation for arbitrary d is not known. This compHcates 
the calculations which involve the wave function. However, since in the con- 
figuration space all divergences are located at r = 0, it is possible to extract 
divergent contributions in the form of 4''^{0)/e. Thus, the unknown overall 
factor ip'^{0) will be present in the intermediate results. In the final result the 
singularities 1/e cancel and we can replace ip{0) by its 3-dimensional value. 

In rare cases we need, in addition to the value of the configuration space 
wave function at the origin, also the behavior of its first derivative. This, 
however, is related to tp{0) by the Schrodinger equation. To demonstrate this, 
we consider the Schrodinger equation in the configuration space. Since for S 
states ijj{r) has no angular dependence, we have 

-— + V(r)] ^P{r) ^ -— (dl + '^—^dr - mVir)] Mr) = 0, (9) 
m J m \ r J 

where V{r) is the Coulomb potential in d dimensions, 

Vir,^-^^. (10, 

We assume that the wave function near the origin behaves like ■(/'(r) ~ ■!/)(0)(l — 
cr''), and find c and P by substituting this form into Eq. (^: 

V^(r ^ 0) V(0) fl - """^ ~ r'+A . (11) 
^ -^W ^ 27r5-'(l + 2e) / ^ ' 

We now would like to demonstrate how divergent integrals are treated in 
practice, using two examples from our hyperfine splitting calculations. First, 
we consider, a logarithmic divergence arising from tree-level Coulomb photon 
exchange, E£l 



In Eq. (12) the matrix element is to be calculated over the ground state wave 



function in d dimensions: 

Although the integrand does not look complicated, the difficulty is that the 
exact form of the wave function if^ir) in d dimensions is not known. We use the 
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fact that it satisfies the d-dimensional Schrodinger equation, Eq. (|^). Using 
that equation we rewrite the integral in Eq. (|l2|) as 

'(P'q)(qp)\ / (47rQ!m)2(p'q)(qp) 



/pp' \(p2-m£;)(p-k)2q2(p'2_„^S)(p'_k')Vkk'' 

' (13) 

where the integration over p, p', as weh as k, k', in the last expression is 
understood. The integral over p and p' receives a divergent contribution only 
from the region where p and p' become simultaneously infinite. Therefore, a 
single subtraction is sufficient to make this integral finite. It is convenient to 
subtract from (|l3|) the following expression: 

(47ram)2(p'q)(qp) 



(p2 — mi?)2q2(p'2 _ mE) 



(14) 



After the subtraction is done, two nice features emerge. In Eq. ( |l4| ) the inte- 
gration over k, k' factorizes and leads to ip'^ (0) times a two- loop integral, which 
can be easily calculated for arbitrary d. On the other hand, the difference be- 
tween the last integral in Eq. (|l3|) and the integral in Eq. ([l^ ) is finite and can 
be calculated for d — 3 using the explicit form of the wave function, Eq. (0) . 

We note that the counterterm (ij) is constructed in such a way that the 
difference mentioned above vanishes for the ground state. This can be easily 
seen by integrating over k, k' in E q. (1^ and noticing that the p, p'-dependent 



terms in the denominator of Eq. (|l4) coincide (up to a normalization factor) 
with the three-dimensional ground state wave functions in the momentum 
representation. 

Finally, for the Coulomb exchange contribution one finds 

^cEms = ^i^'iO) - 41n(ma) - , (15) 

Our second example is a linearly divergent integral arising from the tree- 
level magnetic photon exchange, E3 



Shifting the integration variable p ^ p + k we find that the p-integral in the 
last term is scale-less. In dimensional regularization such integrals vanish. The 
first term in Eq. ( [l6| ) is finite in three dimensions. We obtain 

(p^) ^mE^^iO). (17) 
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3.3 Hard scale contributions 



Perhaps the main advantage of using dimensional regularization in bound state 
calculations is the possibility of computing the hard corrections analytically. 

The contributions of hard scales arise from virtual momenta of order of 
the electron mass. They can be calculated by considering the on-shell e+e^ 
scattering amplitude exactly at the threshold, i.e. for zero relative velocity of 
the incoming electron and positron. Since the amplitude is independent of the 
incoming spatial momenta, it gives rise to four-fermion operators in the low- 
scale Lagrangian or, equivalently, to the S{r) terms in the effective quantum 
mechanical Hamiltonian. 

Technically, this calculation is similar to the derivation of the matching co- 
efficient pOJj-the vector quark-antiquark current in QCD and its NRQCD coun- 
terpart, yfes 

An arbitrary Feynman integral which contributes to the hard scale part of 
the calculation can be written as 



d^fci d^k2 1 



/(ai, ...ag) I /r,_\n /n—\D aai qa2 qas qa4 qa^ qa^ qaj qas qag ^ (1^) 



where 



Si^kl S2^kl Sz = {k^-k2f, Si = kl + 2pk^, 

S5 = kl + 2pk2, Se ^ kf ~ 2pki, 87 ^ kj - 2pk2, 

Ss = (fci - k2 f + 2p{ki - k2), Sg = (fci - k2)^ - 2p(fci - fca), (19) 



and ai, . . . , ag are integers. In practice we find diagrams with only at most 6 
different propagators, so that at least 3 exponents Oi are zero. 

A method for computing these integrals is analogous to the 2-loop virtual 



corrections to quark decays. We have described our algorithm in Section 2.1 



4 Conclusions 

In this talk we have presented technical aspects of our recent calculations of 
two- loop corrections to semileptonic quark decays and positronium properties. 
The methods developed in the context of quark decays, in particular the al- 
gorithm for solving systems of recurrence relations, have found applications 
in positronium physics. We hope that further progress can be made in such 
atomic physics calculations using the tools developed for high energy problems. 
In particular, applying dimensional regularization should facilitate so far in- 
tractable calculations, such as two-loop corrections to positronium decays. 
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